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ABSTRACT 

We measure the large-scale clustering of Mgll AA 2796,2803 absorbers with respect to a population 
of luminous red galaxies (LRGs) at z ~ 0.5. From the cross-correlation measurements between Mgll 
absorbers and LRGs, we calculate the mean bias of the dark matter halos in which the absorbers reside. 
We investigate possible systematic uncertainties in the clustering measurements due to the sample 
selection of LRGs and due to uncertainties in photometric rcdshifts. First, we compare the cross- 
correlation amplitudes determined using a flux-limited LRG sample and a volume-limited one. The 
comparison shows that the relative halo bias of Mg II absorbers using a flux-limited LRG sample can 
be overestimated by as much as ^ 20%. Next, we assess the systematic uncertainty due to photometric 
redshift errors using a mock galaxy catalog with added redshift uncertainties comparable to the data. 
We show that the relative clustering amplitude measured without accounting for photometric redshift 
uncertainties is overestimated by « 10%. After accounting for these two main uncertainties, we find 
a 1-CT anti-correlation between mean halo bias and absorber strength Wr(2796) that translates into a 
1-cr anti-correlation between mean galaxy mass and Wr(2796). The results indicate that a significant 
fraction of the Mgll absorber population of Wr{2796) = 1 — 1.5 A are found in group-size dark 
matter halos of log Mh < 13.4, whereas absorbers of Wr(2796) > 1.5 A are primarily seen in halos 
of log Mh < 12.7. A larger dataset would improve the precision of both the clustering measurements 
and the relationship between equivalent width and halo mass. Finally, the strong clustering of Mg II 
absorbers down to scales of ~ 0.3 Mpc indicates the presence of cool gas inside the virial radii of 
the dark matter halos hosting the LRGs. 

Subject headings: Quasars: absorption lines — Cosmologyitheory — dark matter — galaxies:evolution 
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1. INTRODUCTION 

Characterizing the structure and evolution of the cold 
gas in dark matter halos is a key element in current mod- 
els of galaxy formation (e.g., Keres et al. 2005). The 
fraction of cold and hot gas within dark matter halos 
and the rate at which gas is being accreted are essential 
to our understan ding of disc and star formation (e.g., 
iDekel et"al]l2009f ). Extended gaseous envelopes around 
galax ies were first predicted several decades ago (jSpitzed 
[19561) . Observations of HI maps around local galaxies 
(;e.g.. lThilker"eran[2?)nirDovle et al.l[2?)n5h and compar- 
isons of galaxies and QSO a b sorption-line systern s (e.g. 
Bergeron fc Stasihskal Il986t iLanzetta fc BowenI Il990l: 



Steidel et all 119941 : iChen et all 120011 : iChen fc Tinker! 
20081 ) have indeed shown the presence of extended cool 
gas (T - lO'' K) out to 50 - 100 h'^ kpc radii. The 
physical mechanism that explains the origin of the ex- 
tended cold halo gas is, however, unclear. Some of the 
most common scenarios are (i) outflows from starburst 
systems (e.g.. iBond et al.l l 2001[): Cii) strip ping from the 
accretion of gas-rich satellites (|Wangj[T993() ; (ni) cold gas 
bound to substructur e within the host dark halo (e.g., 
ISternberg et all [2002( 1. and (iv) a two-phase medium 
composed of cold and hot gas (|Mo fc Miralda-Escudd 
[19961: lAMkr fc Bullockl[200l . 

A potential probe of the cold halo gas is the Mgll 
AA 2796,2803 absorption features commonly seen in 
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the spectra of background QSOs. These absorbers are 
thought to originate in photo-ionized gas of tempera- 
ture T ~ 10* K and to trace high-column density HI 
clouds of neutral hydrogen column density A^(HI) 
IQis 1022 pj^-2 (jBergeron fc Stasihska|[T986l : iRao et all 
l2006( l. This large associated HI column density indi- 
cates that Mgll a bsorbers arise in halo gas around in- 
dividual galaxies (jDovle et al.ll2005[ ). This is also sup- 
ported by the presence of luminous galaxies at pro- 
jected distances p = 50 — 100 kpc from known 
Mgll absorbers (iBerger onI Il986[ ILanzetta fc BowenI 
19901: iBergeron fc Boissc 19 9lt ILanzetta fc Bowenlll992l: 



Stcidcl & Sargent 1992: Stcidcl ct al."l994l 'ZibettLeialJ 
2007: Nestor ct al. 2007; Kacprza k ct al ,, 2008{k 

In addition to the classical gas accretion scenario for 
the origin of Mg II absorbers at larger galactic radii, there 
is a competing scenario that has gained substantial atten- 
tion recently. In this new picture, strong Mgll absorbers 
of rest-frame absorption equivalent width Wr(2796) > 1 
A ori g inate in starburst driven out fl ows (e.g..lBoiid et al.l 
[200l iMenard fc Chelouchd [20081 : iWeiner et all I2008D . 
Under this scenario, the Mgll absorbing gas orginates 
in the cold outflowing material surrounding starburst 
galaxies. An interesting recent finding is a strong correla- 
tion between dust ex tinction E{B — V) and Wr(2796) by 
IMenard et all (|2008[ ). While the observed E{B - V) vs. 
Wr(2796) is consistent with the expectation of the star- 
burst scenario, this observation is also expected if the 
Mgll absorbing galaxies exhibit a meta llicity gradient 
commonly seen in regular galaxies (e.g., IZaritskv et al.l 
van Zee et al.l [l"998[ l. Dense clumps in starburst 
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driven outflows are expected to contribute to some frac- 
tion of the observed Mgll absorbers, but the significance 
of this fraction and how the fraction varies with Wr (2796) 
are both uncertain. 

As a first step toward a quantitative understanding 
of the physical origin of the Mgll absorber population, 
we are carrying out a cross-correlation analysis of Mg II 
absorbers with photometrically identified luminous red 
galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS; 
lYork et al.l[2000l) . The primary goals are (1) to deter- 
mine the clustering amplitude of Mgll absorbers and 
(2) to examine how the clustering amplitude depends on 
absorber strength 1^^(2796). The clustering amplitude 
of Mgll absorbers is determined based on their cross- 
correlation signals with LRGs on projected co-moving 
distance scales of = l—30h~^ Mpc. Because the mean 
halo"^ mass of the L RGs can be calcula t ed from their clus- 
tering signal (e.g . IZheng e t aLll2008t iBlake et all l2008l : 
IWake et all 120081 : iPadman abhan et al.ll2008f l the cross- 
correlation amplitude of LRGs and Mg II absorbers pro- 
vides a statistical estimate of the mean mass of the host 
dark matter halos. A similar study was published by 
Bouche et al. (2004, 2006), in which the authors at- 
tempted to constrain the mean halo mass of Mgll ab- 
sorbers using a flux-limited (i' < 21 mag) sample of 
LRGs found in SDSS (see also Lundgren et al. 2009). 

Our analysis differs from others in two important as- 
pects. First, we measure the clustering amplitude of 
Mgll absorbers using a volume-limited (instead of flux- 
limited) LRG sample. A flux-limited selection criterion 
forms an inhomogeneous sample of LRGs, excluding pro- 
gressively more intrinsically fainter LRGs (and hence 
lower-mass halos) toward higher redshifts. Such inho- 
mogeneous samples of LRGs over a broad redshift range 
(z = 0.35 — 0.8) include an inherent uncertainty in the 
estimated mean halo mass that is difficult to assess, but 
this systematic bias has not been addressed in previous 
studies. We have directly compared the clustering ampli- 
tudes between using a fiux-limited and a volume-limited 
LRG sample. We will show in the following sections that 
the relative halo bias of Mg II absorbers may have been 
overestimated by as much as ~ 20% in previous studies. 

Second, the LRGs in SDSS have been identified us- 
ing photometric redshifts (zphot) that have associated 
uncertainties relative to spectroscopic redshifts (zspec) 
of gp = Izphot - ZsD ec l/(l + ^sDoc) « 0.03 at i' « 19 
()Collister et al.l l2007l : lOvaizu et all l2008f l. At fainter 
magnitudes, (Tq increases steeply to ctq > 0.1 at i' > 20.7 
(CoUister et al. 2007). It is clear from these studies that 
the redshift uncertainties of the LRG sample vary from 
Cz = <To{l + z) = 0.04 for galaxies of i' < 19 at z = 0.35 
to = 0.18 for galaxies of i' w 21 at z = 0.8. On the 
other hand, Mgll absorbers are identified in QSO spec- 
tra with a redshift precision better than ct^ w 0.0004, 
corresponding to roughly half of the width of a resolu- 
tion element. The dependence of photometric redshift 
errors on galaxy brightness and redshift are expected to 
introduce additional systematic uncertainties in the esti- 
mate of the Mgll-LRG cross-correlation amplitude. To 
reduce the systematic bias due to photometric redshift 
errors, we first restrict our analysis to including only 

^ We define halo as a region of overdensity 200 with respect to 
the mean mass density in the universe. 



galaxies brighter than i' ~ 20**. This allows us to main- 
tain an LRG sample of higher redshift precisions. In ad- 
dition, we assess the systematic uncertainty using a mock 
galaxy catalog with and without an added ctq x (1 + z) 
redshift perturbation. We will show in the following sec- 
tion that relative clustering amplitude of Mg II absorbers 
measured without accounting for photometric redshift 
uncertainties may have been overestimated by w 10%. 

This paper is organized as follows. We first describe 
our data samples, Mgll and LRGs in Section 2. In 
Section 3.1, we describe the method used to calculate 
the two-point correlation function and associated er- 
rors. Cross-correalation results are presented in section 
3.2, including a discussion of the effects of photomet- 
ric redshift uncertainties in the calculation. In Section 
4, we summarize the routines adopted for calculating 
the bias and mean halo mass of absorbers. Finally, 
we discuss the results of our analysis in Section 5. A 
more extensive halo occupation distribution analysis of 
the Mgll absorbers will be the subject of a forthcom- 
ing paper. We adopt a ACDM cosmology, f^M = 0.25 
and = 0.75, with a dimensionless Hubble parameter 
h = 77o/(100 km s"^ Mpc"^) throughout the paper. All 
masses are expressed in units of Mq unless otherwise 
stated. 

2. DATA 

2.1. Mgll Catalog 

The Mgll absorbers catal og is based on an ex tension 
of the SDSS DR3 sample (jProchter et al.l I2006D to in- 
clude DR5 QSO spectra (J.X. Prochaska, private com- 
munication). This sample of Mgll absorbers has a 
95% completeness for absorbers of VF,.(2796) > 1 A. 
The catalog contains 11,254 Mgll absorbers detected at 
ZMgii = 0.37 - 2.3 along 9,774 QSO sightfines. From 
this sample, we selected absorbers with a separation of 
at least 10,000 km/s from the QSO redshift to avoid con- 
tamination by associated absorbers of the QSOs. We 
excluded absorbers that fall outside of our survey mask 
(see Section 3.1.1) and we limited ourselves to Mgll dou- 
blets in the redshift range of interest, z = 0.40 — 0.70, 
defined by the LRG sample. Note that none of the ab- 
sorbers listed in the Prochter et al. catalog are found in 
the Lya forest of their respective quasar. The procedure 
described above yielded 1,158 absorbers of Wr.(2796) > 1 
A. This contitutes our primary Mgll catalog. The short 
dashed histogram in Figure [1] shows the redshift distri- 
bution of our primary Mgll sample. 

2.2. LRG Catalog 

LRGs are tracers of the large-scale structures in the 
universe. Their clustering properties are well-known, 
and t hey reside in mas s ive halos of Mt > 10^'^ /t~^ M(^ 
(e.g. IBlake et a l.' '2008; Zheng elalj l2008t IWake et alJ 
120081 : IPadmanabhan et al. 2008 ) . Because they are lumi- 
nous and have strong 4,000-A break, their photometric 
redshift can be more reliably estimated than blue star- 
forming galaxies. 

We used the Mega Z-LRG catalog (hereafter MegaZ) 
(iGollister et aLlIMm as our initial LRG sample. MegaZ 

* The volume-hmited selection criterion is applied after adopting 
this brightness limit. 
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is a photometric redshift catalog of approximately 10^ 
LRGs found in the SDSS DR4 imaging footprint. The 
catalog covers more than 5,000 deg^ in the redshift 
range 0.4 < z < 0.7. The LRGs are selected follow- 
ing a serie s of cuts in a multidimensional color dia- 
gram ^e.g.. IScranton et all l2003l : lEisenstein et all l2005l : 
ICollister et al.ll2007[ ). The photometric redshifts of LRGs 
are determined using artificial neural networks (ANNz) 
and the 2SLAQ traini ng set of 13,000 ava ilable LRGs 
spectroscopic redshifts (jCannon et al.|[2b06[ ). 

We applie d a set of addition al color selection criteria 
suggested bv lBlake et al.l (|2007( ). These modified criteria 
yielded less contamination by blue galaxies a nd were used 
in clu stering and halo occupation analyses (jBlake et al.l 
I2008D . We further limited ourselves to galaxies with i' < 
20 for a higher photometric redshift precision. We de- 
fined three LRG samples for our analysis. The first was a 
flux-hmited sample of i' < 20 LRGs at z = 0.4-0.7. This 
primary sample covers the entire redshift range offered by 
the initial LRG sample. The second was a volume-limited 
sample of M,> - 5 log h < -22 at z = 0.45 - 0.60. The 
redshift range was selected to provide the largest num- 
ber of LRGs available under a uniform minimum rest- 
frame absolute magnitude selection criterion. Extending 
to lower or higher redshifts would result in a significant 
reduction in the sample size. The third was a flux-limited 
sample of i' < 20 at z = 0.45 — 0.60 to be directly com- 
pared with the volume-limited subsamplc from the same 
redshift range. 

Our sample definition was motivated by the knowl- 
edge that a flux-limited sample contains galaxies that 
are progressively fainter at lower redshifts, whereas a 
volume-limited criterion identifies a uniform sample of 
galaxies that occupy the same luminosity interval at 
different redshifts. In addition, the cross- and auto- 
correlation calculations are evaluated at different redshift 
in a flux-limited sample. The redshift number density of 
LRGs peaks at z '-^ 0.45, while the Mgll redshift dis- 
tribution is flat. This implies that the mean redshifts 
of the cross-correlation and auto-correlation calculations 
are different, introducing additional uncertainties in the 
estimated mean halo bias. This problem is alleviated 
in a volume-limited sample because the cross- and auto- 
correlation terms have similar redshift evolution (see Fig- 
ure 1). Comparing the correlation functions determined 
using different subsamplcs allowed us to evaluate possible 
systematic uncertainties due to sample selections. 

The limiting absolute magnitude of the volume-limited 
sample was determined by calculating the absolute mag- 
nitude of our faintest galaxies (i' = 20.0) at the limiting 
redshift z = 0.60 

20.0 - A/, = DM(z) + fc(z) , (1) 

where DM(z) is the distance modulus and fc(z) represents 
the fc-correction. At z = 0.60, the limiting magnitude is 
Mi — 5 \ogh = —22 and we found roughly 197 K galaxies 
more luminous than this limit. A quantitative descrip- 
tion of the three LRG samples including the number of 
objects and redshift range can be found in Table [1] Fig- 
ure [1] shows the redshift distributions of the three LRG 
samples and the top panel of Figure [2] shows the magni- 
tude distributions of our flux-limited LRG sample and of 
the initial LRG sample. 
It is notoriously difficult to determine the photometric 
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Fig. 1. — Redshift distributions of the initial McgaZ LRG sample 
and of the three LRG samples used in our correlation studies. The 
initial MegaZ catalog is shown in the dotted line. The solid line 
shows the flux-limited sample (z = 0.40 — 0.70) and the dashed 
line shows the volume-limited sample (z = 0.45 — 0.6). The redshift 
distribution of Mg II absorbers in our primary sample is overplotted 
in short-dashed line for comparison. 

redshift of faint galaxies. The bottom panel of Figure 
[5] shows a subsample of the photometric redshift errors 
taken for 1,000 overlapping galaxies between MegaZ and 
the SDSS photometric redshift table. The photoz's er- 
rors shown in this figure are taken in the "photozcc2" 
table found on the SDSS skyserver archive (jOvaizu et al.l 
I2008| ). The photometric redshift errors ob tained by the 
MegaZ team are consistent with the ones in lOvaizu et al.l 
()2008( ) to within measurement uncertainties. However, 
it is clear from Figure [2] that for i' > 20, the uncer- 
tainties in the photoz's increase rapidly. For this reason, 
we decided to restrict ourselves to LRG candidates with 
i' < 20. This procedure yielded a total of 962,216 MegaZ 
objects satisfying our additional selection criteria. This 
constitutes our primary LRG catalog. 

3. MGII ABSORBERS-LRG CROSS-CORRELATION 
3.1. Method 

We used the lLandv fc Szalavl (|1993[ ) (LS93) minimum 
variance estimator to calculate the projected two-point 
correlation statistics between Mgll absorbers and LRGs. 
Using the LS93 estimator, the real-space correlation 
function can be calculated following 

, N Di^Dg - DaRg - DgRa + ^a^g 

^) ^ (^) 

where R and D are random points and data, the sub- 
scripts a and g refer to absorbers and galaxies, rp is the 
projected comoving separation between two objects on 
the sky and tt is their distance parallel to the line of 
sight. Note that the above two-point estimator has been 
successfully used in previous correlation study based on 
QSO absorb ers and galaxy data col lected from a smaller 
area survey (jAdelberger et alll2003[ ). In practice, we cal- 
culated the projected two-point correlation statistics by 
summing all pairs along the sightline 

^pirp) J^^{rp,Tr)d'K (3) 
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Fig. 2. — Top panel: Magnitude distribution of the LRGs satis- 
fying our selection criteria 9.6 X 10^ LRGs). The line shows the 
i' magnitude distribution. Only galaxies with i' < 20 are included 
in our LRG sample. Bottom panel: Photometric redshift error ver- 
sus i ' magnitude f o r a su bsample of 1,000 galaxies in SDSS taken 
from lOvaizu et all l)2008I V The photometric errors are listed in the 
SDSS Table "photozerr2" under "CC2". Black diamond points 
with error bars represent the mean and dispersion calculated in 
bins of Ai' = 1 mag. 

and inside of redshift limits I of our three LRG samples. 
Another commonly used estimator, l-|-^(rp,7r) = ^"^"^ 

has also been used in the initial stages of this work for 
comparison and validation purposes. Even though this 
estimator is easier (no random absorbers) and faster (two 
terms to compute instead of four), the variance of Wp is 
larger for this estimator than for LS93. For this reason, 
we employed the LS93 estimator for all Wp calculations. 

We divided the pairs into eight rp bins covering the 
range 0.2 — 35 Mpc. The bins were equally separated 
in logarithm space. The bin size and the rp value of the 
inner most bin were determined in such a way that at 
least 10 DD pairs were found in that bin. The upper 
limit of 35 Mpc was chosen to be a few times smaller 
than the size of the jackknife cells (see section 3.1.3). 

3.1.1. Survey mask 



For the calculation of the two-point statistics (follow- 
ing equation [2|) both data and randoms were distributed 
over the exact same survey mask. Different masks for 
LRGs and Mgll would alter the shape and amplitude of 
the correlation signal in an undesirable fashion, especially 
at large separations. To make sure the sky coverage was 
the same for Mgll, LRGs and their randoms, we used 
the lowest common denominator for all : the DR4 spec- 
troscopic sky survey mask. Indeed, Mgll absorbers were 
taken from the SDSS DR5 spectroscopic sample which in- 
cludes the QSO sightlines inside the DR4 spectroscopic 
sky. The LRGs were found in the DR4 imaging sky which 
encompasses the DR4 spectroscopic sky. The main dis- 
advantage of using this mask was that the number of 
rejected objects falling outside of the mask was large, 
«40%. 

To be able to determine which Mgll absorbers and 
LRGs fell inside this DR4 spectroscopic mask, we used 
the mask catalog pro vided by the NYU Value-Added 
Galaxy Catalog team iBlanton et al.l (j2005[ ) . The mask 
corresponds to the angular selection function describ- 
ing the completeness of the SDSS spectroscopic across 
the sky. It is defined by spherical polygons. The com- 
pleteness quantifies the fraction, inside each polygon, of 
galaxies with spectrocopic redshift. We used the angu- 
lar selection function of the spectroscopic SDSS DR4 sky 
(sdssdr4safe0res6d) . 

3.1.2. Generating random Mgll absorbers and LRGs 

The right ascension and declination of the LRGs 
were randomly selected over the DR4 spectroscopic 
sky using the function ransack available in the 
Mangle software pa ckage (jHamilton fc Tegmarkl l2004l : 
ISwanson et al.l[200a ) . The redshifts of the random galax- 
ies were determined by sampling the redshift distribution 
of the LRG dataset. Determining the sky positions and 
redshifts of the random absorbers was not a straightfor- 
ward process. The angular selection function of quasars 
follows the DR4 spectroscopic mask, but the mask defin- 
ing the positions where absorbers can be found is lim- 
ited to the actual coordinates of the QSOs themselves. 
Thus, the random absorbers must be distributed ran- 
domly among fixed QSO sightlines for which SDSS spec- 
tra are available. Assigning random Mgll this way elim- 
inates any undesired bias due to the intrinsic clustering 
of QSO sightlines. We id entified these sightlines from 
the lSchneider et all ([200l SDSS DR5 QSO catalog by 
selecting all QSOs falling inside our DR4 spectroscopic 
mask with redshifts large enough to allow for the de- 
tection of Mgll absorbers inside the redshift range of 
interest for our calculations (z = 0.40 — 0.70). We found 
a total of w 5.5 x 10** sightlines. The (ra,dec) positions 
of the random absorbers were determined by randomly 
selecting the coordinates of these sightlines. Redshifts 
were selected randomly from a top hat probability dis- 
tribution function over z ~ 0.4 — 0.7. 

The number of random LRGs, Rg, and the number 
of random Mgll absorbers, Ra, were determined af- 
ter running convergence tests with increasing number 
of randoms. We varied the number of each random 
set (Mgll and LRGs) independently until the measured 
Wp approached asymptotic values. The numbers of ran- 
doms used in our Wp calculation were Ra = 10'^ and 
i?g = 4 X 10^ for Mgll and LRGs, respectively. It is 
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important to note that increasing the number of random 
Mgll absorbers to a very large number in comparison 
with the number of available QSO sightlines (~ 55K) 
would not lead to a rapid convergence because we would 
simply be overcounting the same D^R^ pairs and no ad- 
ditional information would be gained. 

3.1.3. Relative contribution of cosmic variance and 
photometric uncertainty to Wp errors 

We looked at two independent sources contributing to 
the error bars on Wp : cosmic variance and photomet- 
ric redshift errors. Photometric redshift errors are ex- 
pected to affect the measurements of Wp^s in two differ- 
ent ways. First, redshift uncertainties are expected to 



increase the noise in the 



measurements due to un- 



certainties in the object positions. Second, redshift un- 
certainties are expected to systematically alter the Wp 
measurements to lower values due to an inherent sam- 
ple selection bias. A galaxy sample selected based on 
imprecise photometric redshifts contain galaxies from a 
broader redshift range than one selected based on pre- 
cise spectroscopic redshifts, and consequently reducing 
the observed correlation signal by including additional 
uncorrelated pairs. In this section, we address random 
errors in the Wp measurements due to cosmic variance 
and object distance uncertainties. We defer the discus- 
sion on the systematic errors of Wp due to the sample 
selection bias to § 3.2.2. 

We estimated the cosmic variance using the jackknife 
resampling technique. The sky was separated into N = 
192 (see section 3.2.1 for a justification of this number) 
cells of roughly equal survey area. The cosmic variance 
for each point rp ^ corresponds to the ith-diagonal el- 
ement of the covariance matrix was calculated using a 
jackknife resampling technique : 



COVK,w. 



N 



1^ 



fc=i 



{w^ ~uh){w'; ~wj) (4) 



and k represents the iteration in which box k was re- 
moved. The mean Wi was calculated for bin i over all 

The impact of the large uncertainties due to pho- 
tometric redshifts on the size of the Wp errors is not 
taken into account by the jackknife resampling tech- 
nique alone. Pr evious works (e.g., Bouclie ct al. 2006t 
iBlake et al.ll2008[ ) have focused primarily on cosmic vari- 
ance in the calculation of the errors bars, but have not ad- 
dressed additional random noise due to photometric red- 
shift errors which are quite large (ctz ~ 0.05 for i' < 20) 
compared to the redshift (z = 0.5). 

To account for the independent contribution of pho- 
toz's uncertainties on the final Wp error bars, we gen- 
erated 100 independent realizations of the MegaZ cata- 
log. For each one of them, we resampled the redshift of 
each individual galaxy according to a normal distribution 
N(zpiiot, o'z), where Zphot is the photometric redshift and 
az is the photometric redshift error of each galaxy. In 
this cas e, we followed the erro r function a, = 0.03(1 + z) 
found in lColhster eTall (|2007t ) to assign photometric red- 
shift errors to galaxies. A new mock LRG sample was 
then established according to the criteria discussed in § 
2.2. We calculated Wp for each one of these realizations 
and assigned the error contribution from redshift uncer- 



tainties to the final error budget of Wp by calculating the 
dispersion among these 100 independent realizations. 

We found that the size of the Wp error bars was dom- 
inated by cosmic variance. The contribution of photo- 
metric redshift errors to the Wp error budget was small, 
at most 20% for the smallest separation bin. It is neg- 
ligible (~ 1%) at large separations where we calculated 
the relative clustering strength and absolute bias. For 
this reason, we decided to adopt the cosmic variance <7i 
as the error on Wp and neglect the contribution of the 
photometric redshifts. 

3.2. Results 

This section addresses the cross- and auto-correlation 
results for the three LRG samples and the effect of pho- 
tometric redshifts on the clustering amplitude. For each 
one of our three LRG samples, we considered three sub- 
samples of Mgll absorbers: weak Wr(2796) = 1 — 1.5 A, 
strong Wr{2796) = 1.5 — 5 A, and all absorbers. A de- 
scription of each correlation calculation can be found in 
Table [T] We also include the number of DD pairs (equa- 
tion [J) in the first bin and the number of data found in 
each subsample in columns (6), (7), and (8). 

3.2.1. Cross- and auto- correlation results 

Figure [3] shows the cross- and auto-clustering results 
for the flux- and volume- limited samples of LRGs. For 
each one of these samples, the results are shown for the 
three Mgll subsamples in gray. We overplot the LRGs 
auto-correlation results in black. A clear feature of Fig- 
ure [3] is that, for the three LRG samples considered, the 
clustering amplitude of weak absorbers is systematically 
higher than for strong ones. It is also interesting to note 
that after accounting for systematic bias due to redshift 
uncertainties (see § 3.2.2), weak Mgll absorbers appear 
to be unbiased (sharing the same clustering amplitude) 
with respect to LRGs in both flux-limited and volume- 
limited samples. 

For each auto- and cross-correlation calculations, we 
estimated the correlated uncertainties between different 
Tp's using the covariance matrix (see equation[5]) and the 
normalized correlation matrix 



COY {w„Wj) 



(5) 



We performed a convergence test on the number of jack- 
knife boxes. We made sure that the off-diagonal elements 
of the correlation matrix varied by less than 10% after 
doubling the number of boxes. We plot in Figure [4] the 
normalized correlation matrix p['Wi,'Wj) for all calcula- 
tions, including the LRGs auto-correlation shown in the 
fourth column. Adjacent bins are strongly correlated for 
the auto-correlation functions of LRGs and, in all cases 
of the LRG-Mgll cross-correlation measurements, bins 
with large Vp are more correlated. 

3.2.2. Effect of photometric redshifts on the amplitude of 

clustering 

An important effect of photometric redshift errors is 
the sample selection uncertainty. Photometric redshift 
uncertainties not only affect the precision of Wp, but also 
lower its accuracy that one wo uld normally get fr om us- 
ing spectroscopic redshifts (e.g.. lBrown et ari2008D . This 
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TABLE 1 

Summary of the cross- and auto-correlation calculations 









Wr.n.in{2796)" Wr,max(2796) 


Number 


Number 


Number 


Sample 


^min 


■2max 


^ (A) A 


DD pairs(istbin)'' 


LRGs 


Mgll 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


Volume-limited sample 


V-weak(VW) 


0.45 


0.60 


1.0 1.5 


21 


197,968 


279 


V-strong(VS) 


0.45 


0.60 


1.5 5.0 


13 


197,968 


257 


V-all(VA) 


0.45 


0.60 


1.0 5.0 


34 


197,968 


536 


Flux-limited samples 



Fl-weak(FlW) 


0.45 


0.60 


1.0 


1.5 


47 


517,549 


279 


Fl-strong{Fls) 


0.45 


0.60 


1.5 


5.0 


35 


517,549 


257 


Fl-all(FlA) 


0.45 


0.60 


1.0 


5.0 


82 


517,549 


536 


F2-weak(F2W) 


0.40 


0.70 


1.0 


1.5 


88 


618,086 


541 


F2-strong{F2S) 


0.40 


0.70 


1.5 


5.0 


68 


618,086 


617 


F2-all(F2A) 


0.40 


0.70 


1.0 


5.0 


156 


618,086 


1158 



LRGs auto-correlation 



V-LRGs(VG) 0.45 0.60 - - 6,977 197,968 

Fl-LRGs(FlG) 0.45 0.60 - - 43,093 517,549 

F2-LRGs(F2G) 0.40 0.70 - - 55,062 618,086 



^ Wr is the rcst-framc equivalent width. 

^ The first hin is centered on 0.312 h^^ Mpe 

needs to be accounted for in the bias calculation. One 
can think of photometric redshift errors as adding uncor- 
related galaxies in the calculation, and the overall effect 
is an unwanted widening of the galaxy redshift distri- 
bution. This effect results in a systematic error in the 
clustering signal and the amplitude of Wp is compara- 
tively smaller than expected in calculations using only 
spectroscopic redshifts. 

To estimate the attenuation of the clustering signal 
due to photometric errors, we performed several tests on 
mock LRG galaxy distributions. Our mock catalog is 
produced by populating the halos in an N-body simula- 
tion with a halo occupation f unction deterrnined from the 
spectroscopic LRG sample in lZheng et al.l ()2008[) (specif- 
ically, the "faint" sample). The simulation is large, 1 
Gpc/h on a side , and represents ou r fiducial cosmology 
at z - 0.5 (see iTinker et all l2008l for more details on 
this simulation). Since the box size is not big enough to 
cover the redshift range of our study, we mirror-imaged 
three identical copies of the box along the redshift direc- 
tion. The mock LRG catalog has precise positions and 
therefore mimic a spectroscopic sample. 

Three tests were performed. First, we auto-correlated 
the mock LRG sample without applying any modifica- 
tions to their redshifts. This first case mimics the results 
one would get by auto-correlating a sample of LRGs with 
spectroscopic redshifts (the Whh case in Figure 6). Sec- 
ond, we cross-correlated the mock LRG sample with one 
that involved perturbed redshifts. The perturbed red- 
shifts were generated by sampling a normal distribution 
centered on the input redshift in the mock catalog within 
(Tz = 0.03(l+z). That is the Whh' term in Figure[Sl which 
is in analogous to our cross-correlation measurements be- 
twen absorbers with spectroscopic redshifts and galaxies 
with photometric redshifts. Third, we auto-correlated 
the mock LRG sample with perturbed redshifts {wh'h' in 
Figure 6). This mimics our LRG auto-correlation func- 
tion. 



We then computed the Wp ratios between the differ- 
ent calculations and plotted the results in Figured The 
three sets of symbols correspond to ratios of Wp calcu- 
lated over the same series of bins as in the correlation 
calculations shown in Figure [S] We show, by thick black 
lines, the best-fit amplitude of the three ratios in two 
different regimes : small and large separations. Best- 
fit values are listed in Table [2] with their corresponding 
uncertainties (Ax^ < !)■ Error bars on the ratios them- 
selves correspond to the dispersion obtained among 100 
realizations of a mock catalog with perturbed redshifts. 
The best-fit amplitude (of the light points) at large sep- 
aration enters in the calculation of the relative bias (see 
section 5.1). 

The dark datapoints in Figure [5] show that the cross- 
correlation function between Mgll absorbers and pho- 
tometric redshift identified LRGs is underestimated by 
roughly 20% due to redshift errors. The gray points 
show that the auto-correlation function of these LRGs is 
underestimated by roughly 30%. Finally, the light-gray 
points show that the mean bias calculated based on the 
ratio of absorber-LRG cross-correlation and LRG auto- 
correlation functions is overestimated by roughly 10%. 

Our estimated reduction in clust ering strength of 10% 
is smaller than the value quoted in lBouche et al.l ()2004f ). 
These authors found that, in the case of a gaussian red- 
shift distribution for galaxies, the amplitude of the Mgll- 
LRG cross-correlation is overestimated by 25 ± 10%. 
They used numerical integration and mock galaxy cata- 
log with phometric redshift uncertainties = 0.1 cor- 
responding to size of their redshift interval of interest. 
The discrepancy could be partly attributed to the larger 
photometric redshift uncertainties used by these authors. 

4. DETERMINING THE ABSOLUTE BIAS AND MEAN 
MASS SCALE OF ABSORBER HOSTS 

4.1. Theoretical Framework 

The bias of dark matter halos can be defined as the 
ratio between the clustering of halos (at a fixed mass) 



Mgll & LRGs at z ~ 0.5 



7 



Flux-limited samples 



1000 
100 
10 
lOOd) 



Q- 100 



10 r 



1 



-| 1 I I I I ll| 1 1 I I I I ll| 1 TT 

♦ Wi2796)=1-5 



T 1 — I I I 1 1 III 1 — I I I I I III 1 — I i-t 



11 



■ I III 



1 1 1 11 1 — I — I 1 1 1 1 1 



f 



it 



t 



! 



I Wr(2796)=1-1.5 



♦ W,(2796)=1.5-5 



T ♦ 



M i 




▼ 



<> :: 



-| 1 I I I I III 1 1 I I I I III 1 TT 



'I 



o 
1^ 



o 



N 



O 
CD 

O 

uS 

d 

II 

N 



Volume-limited sample 



1000 
100 
10 
1 



1 1 1 1| 1 — I — I I I 1 1 1| 1 — r 

Wr(2796)=1-5 J 




_l I t I ^ ■ ■ 



■V I 



"T — I I I 1 1 11 1 — r- 



Wr(2796)=1-1.5 



-I ■ ■ 



u 



-I 



< ► :: 



rn-m 1- 



Wr(2796)=1.5-5 



If{j 



_l 



-I 



0.1 



1 



10 0.1 



1 



10 0.1 



10 



o 

CD 



II 

N 



rp [h"^ Mpc] 



Fig. 3. — Cross- and auto-correlation results for the flux and volume-limited samples of LRGs. Prom top to bottom, flux-limited sample 
z = 0.40 — 0.70, flux-limited z = 0.40 — 0.60, and volume-limited z = 0.40 — 0.60. From left to right, all absorbers, weak absorbers and 
strong absorbers. The rest-frame equivalent-width coverage is listed at the top of each column. In all panels, the gray points represent 
the cross-correlation results (Mg II and LRGs) and the black points are for the LRGs auto-corrclation. In the bottom left panel, wc show 
the halo occupation distribution analysis done on the volume-limited sample of LRGs. The solid line shows the full HOD model and the 
dotted line the one-halo contribution (see section 4.2). Note that both cross- and auto-correlation results were corrected for photometric 
redshifts (see section 3.2.2). Error bars represent the l-tr jackknife errors calculated on N = 192 sky boxes of equal area. Note that the 
last cross-correlation points in two of the strong Mgll absorber sample are negative and not displayed in these panels. The weak absorbers 
are essentially unbiased with respect to the LRGs and their clustering is strong, even down to 0.3 h"^ Mpc. 



and the underlying clustering of dark matter, 



(M,r) 



UM,r) 



(6) 



where ^m('') is the correlation function of the dark mat- 
ter itself. At large scales, linear bias holds and b/^ is 
independent of r. In the translinear regime, r < 10 
Mpc, bh has a scale dependence (with respect to either 
^mif) obtained from linear theory or the true non- linear 



clustering). Although the scale dependence of halo bias 
varies with halo mass, over the mass range probed by 
LRGs and Mgll absorbers, lO^^'^^ Mq, the scale 
dependence is nearly independent of mass and will di- 
vide out in the cross-correlation function (jTinker et al.l 
l2009t) . The auto-correlation function of LRGs can then 
be expressed as 



C,(r) ^ bl{M„r)im{r) = bl{M,)fJr)U{r) (7) 
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Fig. 4. — Gray-scale plot of the correlation matrix (see equation [Sjl showing the degree of covariancc between different separation bins. 
The different matrices shown correspond to the same calculations of Figure |3] for iyr{2796) = 1 — 5 A (left), Wr(2795) = 1 — 1.5 A 
(middle-left), and VFr(2796) = 1.5 — 5 A (middle-right), respectively. The fourth column shows the LRGs auto-correlation results. 



TABLE 2 

Estimated reduction factor of the cross-correlation amplitude due 
TO photometric redshift sampling errors 



Measurement 


Vp < 1 h-^ Mpc 


xVd-O.f 


rp > 1 ft-i Mpc 


xVd-O.f 


WMgll-LRG 


0.77 ±0.03 


1.35 


0.79 ±0.02 


0.26 


""LRG-LRG 


0.66 ±0.04 


0.97 


0.71 ±0.04 


0.12 



where Mg is the bias- weighted mean mass scale (see equa- 
tion [T7| of galaxies, hh and hg are the large-scale linear 
biases of dark matter halos and gala xies, and /(r) is the 
scale -dependent bias term (see, e.g.. [Tinker et al.ll2005L 
[20(11) . The cross-correlation is then 

iga{r) = bg{Mg)ba{Ma)fg{r)fa{r)U{r). (8) 

Thus the relative bias of absorbers to LRGs is the ratio 
of the cross to auto-correlation functions, ie, 



which should be close to a constant at large separations. 
We used the measured Wgg{rp) to obtain &g, which we 
utilized to obtain ba- At z = 0.5 and at masses above 
M ~ 10^^-^ h~^MQ, the bias of dark matter halos in- 
creases monotonically with mass, thus the equivalent 
dark matter halo mass can be obtained from inverting the 
bh(M) form ula. We used the following halo bias function 
({Tinker et al. 2009) for halos defined at an overdensity 
of 200 times the background 



(9) 



6(cr) = l-A- 



B(j- 



CCT- 



(10) 
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TABLE 3 

Best-fit power-law f{x) = ax'' results for the cross and 

AUTO-CORRELATIONS. 



Sample a 6' x^/D.O.F" 



vw 


0.337 ±0.11 


-0.835 


1.2 


0.3 


vs 


0.206 ±0.09 


-0.835 


0.394 


0.13 


VA 


0.275 ±0.075 


-0.835 


1.06 


0.27 


FIW 


0.245 ±0.080 


-0.783 


0.54 


0.14 


FIS 


0.194 ±0.070 


-0.783 


4.60 


1.15 


FIA 


0.206 ±0.056 


-0.783 


3.78 


0.76 


F2W 


0.160 ±0.046 


-0.781 


3.47 


0.87 


F2S 


0.114 ±0.047 


-0.781 


0.776 


0.258 


F2A 


0.136 ±0.035 


-0.781 


3.68 


0.92 



LRGs auto-correlation 



VG 0.356 ±0.013 -0.835 ±0.029 11.25 3.75 

FGl 0.273 ±0.008 -0.783 ±0.026 9.72 3.24 
FG2 0.241 ±0.007 -0.781 ±0.026 6.79 2.26 



^ In the case of the LRGs auto-coirelation calculation, both a and 
b are free parameters. For cross-correlation Mgll-LRGs, b is kept 
fixed and corresponds to the best-fit value obtained for the auto- 
correlation. 

per degree of freedom. For auto-correlation calculations, 
D.O.F— 4 and D.O.F — 5 for cross-correlation. 



TABLE 4 

Relative and absolute bias. Bias-inferred masses. 



Sample 
(1) 



direct 
(2) 



power-law 
(^RA> 
(3) 



direct 

(4) 



power-law direct power-law 
(^RA> (logJ\4)DR (logA4)RA 

(5) (6) (7) 



Relative bias 



Absolute bias 



Mean mass' 



VW 


0.77 


±0.28 


0.85 


±0, 


,25 


1, 


.56 ±0.57 


1.72 ±0.51 


12.8! 


fO.5 
-1.1 


13 0+°-" 


vs 


0.54 


±0.23 


0.52 


±0, 


,21 


1, 


.09 ±0.47 


1.05 ±0.42 


12.0- 


fO.8 


11.9+0.8 


VA 


0.67 


±0.19 


0.69 


±0, 


16 


1, 


.36 ±0.38 


1.41 ±0.33 


12.5; 


fO.4 
-0.8 


12 6+°-"^ 


FIW 


0.80 


±0.26 


0.81 


±0. 


,24 














FIS 


0.75 


±0.25 


0.64 


±0. 


,21 














FIA 


0.74 


±0.19 


0.68 


±0, 


16 














F2W 


0.69 


±0.17 


0.60 


±0. 


15 














F2S 


0.44 


±0.17 


0.43 


±0. 


16 














F2A 


0.56 


±0.13 


0.51 


±0, 


12 















^ The subscript DR represents the results obtained with the mean ratio of the two-halo term data 
points and RA is calculated from the ratio of the best-fit power-law amplitudes of the cross- and 
auto-correlations. 

5 is the absolute bias. 

* Note that we did not quote a lower limit on the halo mass for strong absorbers. Indeed, the bias 
vs. halo mass relationship flattens to bh{M) ~ 0.7 for logM^ < 9. When the lower limit on the 
bias is < 0.7, this gives us no lower limit on the halo mass. 
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Fig. 5. — Ratios of cross-correlation signals obtained from 
the simulation boxes to address the impact of photometric red- 
shift uncertainties on the amplitude of Wp. Solid black : ratios 
of lOp's where one of the two datasets has redshift position per- 
turbed (within its error bars) {wf^y) over tOp's obtained where 
both datasets have spectroscopic redshift (ui^^); dashed gray : 
ratios where both datasets have perturbed redshifts (li'jj'h') over 
two datasets with unperturbed redshifts (uifih); Dotted light gray 
: Wp for perturbed-perturbed datasets (lOf^/^/) over perturbed- 
unperturbed (wf^y). The last set of points gives us an estimate of 
the systematic bias in the observed clustering amplitude of Mgll 
absorbers using photometric redshift identified LRGs. Error bars 
show dispersion among 100 realizations of the perturbed datasets 
{h'). The thick solid lines represent best-fit amplitude obtained 
for the one (rp < l.lh~^ Mpc) and two-halo (rp > l.lh~^ Mpc) 
terms. Best-fit values are listed in table [2] The best-fit amplitude 
obtained for the light gray points is used to correct for the addtional 
systematic error introduced by using photoz's in the calculation of 
the correlation functions. 

where cr is the hnear matter variance on the Lagrangian 
scale (radius of the halo in the initial mass distribu- 
tion when 6{p) = {p ~ p)/p ~ 0) of the halo. R = 
(3M/47rp)^/^, and a, b, c, A, B, C are constants (a = 
0.132, b = 1.5, c = 2.4, A = 1.04, B = 0.4, C = 0.99). 

4.2. The bias of LRGs 

We obtained the bias of LRGs through halo occupa- 
tion modeling of the Wgg data and the number density of 
galaxies in the sample. Note that Wgg was corrected for 
the systematic error due to photometric redshift uncer- 
tainties (see § 3. 2.2). Our modeling was similar to the 
one performed in lZheng et al.l ()2008D . which is ba sed on 
the an alytic halo occupa t ion ni odel developed in IZhengj 
(|200l and [Tinker eFall pOOl . The best-fit model is 
shown in the bottom left panel of Figure [31 which yielded 
a of 9.8 using the full covariance matrix. Halo occu- 
pation models separate pairs from galaxies located inside 
the same dark matter halo (one-halo term) and pairs 
from galaxies located in two different halos (two-halo 
term). The one-halo contribution to the clustering ampli- 
tude of LRGs is shown in Figure [31 At small separations, 
{vp < 1 h^^ Mpc) the one-halo term dominates but the 
two-halo term shapes the clustering signal for > 1 
Mpc. Because the analytic model fully incorporates the 
scale-dependent bias of dark matter halos, we can ob- 
tain the linear bias directly from these data. The bias of 
LRGs in our sample is bg = 2.023 ± 0.006. The high pre- 



cision is due to the large volume of the sample, yielding 
an excellent measurement of the amplitude of Wgg in the 
two-halo regime. 

4.3. The relative bias of absorbers 

We calculated the relative bias using equation ([9]) 
on large scales only (> 1.0 Mpc). Two meth- 
ods were employed. In the first case, we fitted the 
cross- and auto-correlation results using a power-law 
model and estimated the relative bias using the ratio 
of the best-fit amplitudes. This is a standard procedure 
that has been commo nly done in previous works (e.g., 
[Davis fc Feeble j 119831 ) ■ However, the powe r-law model 
does not have a physica l justification (e.g., [Blake et al.l 
[20081 : [Zehavi et al.l[200 ^ . It simply provides an adequate 
fit to the data. In the second case, we directly calculated 
the relative bias by taking a weighted mean ratio of all 
points at fp > 1 Mpc. Because the measurements 
and measurement errors vary significantly between data 
points at different r^'s, we adopted the weights w^'s that 
were designed to maximize the significance of the mean 
relative bias (6), 



1=4 



where 



(11) 



(12) 



the index i denotes the rp bin and ai is the associated er- 
ror of Wag,i/wgg^i computed using the error propagation 
technique. 

The best-fit power-law parameters can be found in Ta- 
ble [3l We first determined the best-fit parameters of 
the LRG auto-correlation function by minimizing the 
function that accounts for the correlated errors between 
adjacent bins: 

^ {w -wfCOY-\w -w) . (13) 

w is the model and w is the data vector. Next, 
we adopted the best-fit slope of the LRG auto- 
correlation function for all corresponding Mgll-LRGs 
cross-correlation calculations. For the LRGs auto- 
correlation, the errors on the parameters were deter- 
mined from all values within Ax^ < 2.3 (two parame- 
ters to fit) from the minimum. The relative and absolute 
biases derived from the ratio of the best-fit power-law 
amplitude are denoted by the subscript "RA" in Table 
4. In the cross-correlation cases, the error on the best- 
fit amplitude a corresponds to all models with Ax^ < 1 
from the minimum value. 

For the relative bias derived from the direct ratio " DR" 
of datapoints in the two-halo regime (equation lll[) . we 
excluded all negative datapoints from the calculation. 
The final error was derived using the standard error prop- 
agation technique and we kept the two most dominant 
terms of the expansion. These two terms are at least two 
orders of magnitude larger than any other term in the 
expansion, 

„2 \ " , .2 1 „2 , 

^ 2^ w,- — cr,„ ; -I- 



2^^tJiWjC0V(Wa 



J>1 



Waga'U'ag.j 



(14) 
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Fig. 6. — Panel (a) shows the relative bias of Mgll for two of the three LRG samples. Solid lines wit h triangles are for th e volume- 
limited sample, solid lines with squares for the flux-limited z = 0.40 — 0.70, and open circles are taken from lBouche et al.l II2006I ) for direct 
comparison with the flux-limited sample. For each LRG sample, the points are derived from the direct ratio of the points for rp > 1 h~^ 
Mpc. The power-law technique yielded very similar results. We shifted the results for the flux-limited samples by {b) + 1.5 for more 
clarity. Panel (b) : absolute bias for the volum e -limit ed sample. Panel (c) bias-inverted halo mass derived from the absolute bias. The 
bias-weighted halo mass results of lBouche etaLl (|2006|) are shown in open circles. Error bars along the x-axis represent the binning used. 



where the sum runs over the bins i,j with Vp > 1.1 
Mpc. Wag,i represents the ith bin of the cross-correlation 
function and Wggj is the jth bin of the LRGs auto- 
correlation function. The relative bias values, calculated 
using both methods, arc listed in columns (2) and (3) of 
Tabled] Wc corrected the relative bias values for photo- 
metric redshift errors by using the large scale correction 
factor 0.90±0.02. 



4.4. The absolute bias and mass scale of absorber hosts 

Since the flux-limited sample of galaxies is not based on 
a homogeneous sample of galaxies and would render the 
halo occupation analysis much more challenging, we lim- 
ited our absolute bias calculation to the volume-limited 
data. The absolute biases obtained, for the direct-ratio 
case, are b = 1.56 ± 0.57 for the weak Mgll sample and 
b = 1.09 ± 0.47 for the strong one. Biases derived from 
the power-law technique arc found in column (5) of Table 

H 

We determined the corresponding halo mass in two 
different ways. In the first case, we inverted equation 
pUj) to obtain the halo mass for these absorbers. These 
masses are denoted by log Mh, which we refer to as the 
bias-inverted mass. Using the direct-ratio evaluated ac- 
cording to equation pT|) . we derived logM^ 



io n-fO.4 
^■^■^-0.6 



for the weak absorbers and log Mh = 12.0_e o for the 
strong ones. The lower error bar we quote for the weak 
absorbers is arbitrary. All halos with logM^ < 12 arc 



consistent with our results. This is because the lower 
error bar on the bias was less than 0.7 giving us no con- 
straint on the halo mass. Indeed, b{M) has a minimum 
value of ^ 0.7 and becomes nearly independent of mass 
at log Mil < 9. The corresponding halo masses using the 
power-law method are in column (7) of Table [H In the 
second case, we first solved for the minimum halo mass 
(Mmin) using 



dlogMh 



Mrr. 



dnjMh) 
dlogMh 



b{Mu) (15) 



where (6) is the absorbers bias, dn{Mh) / dlogMh is the 
halo mass function and 



nh 



dlogMh 



M„ 



dnjMh) 
dlogMh 



We then estimated the mass following 



{logMh)^n^^ 



Mrr. 



dlogMh^^logMh 
dlogMh 



(16) 



(17) 



We refer to the mass found using equation p?)) as the 
bias-weighted halo mass. Mmtn corresponds to the min- 
imum halo mass above which there is a single absorber 
per halo. Of course, absorbers are distributed over a wide 
range of masses and the covering fraction is likely to vary 
with halo mass. This method only gives an approximate 
answer. As for the the bias-inverted mass, it is a char- 
acteristic mass obtained by inverting the b{M) relation. 
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Even though these two methods involve different assump- 
tions, they give similar (within 0.1 dex) results for both 
LRG and absorbers masses. We used the bias-inverted 
masses as our LRG and absorbers mass estimate. 

The relative and absolute biases, calculated with the 
direct ratio technique, and their corresponding halo 
masses are shown in Figure IHl Panel (a) shows the 
relative bias obtained for the volume- and flux-limited 
(z = 0.40 — 0.70 'I samples and allow s for a direct com- 
parison with the iBouche et al.) ()2006D results; (b) shows 
the absolute bias derived from the halo occupation anal- 
ysis. The bias -inverted halo masse s are shown in panel 
(c) along with iBouche et all (|2006( ) bias-weighted mass 
estimates. We did not quote a lower limit for the halo 
mass when the lower limit on the absolute bias is < 0.7. 

We obtained similar relative biases for the flux-limited 
sample a.t z = 0.40 — 0.70 as Bouche et al., who cor- 
rected their relative bias by 20% to account for photot- 
metric redshifts. This large correction is expected since 
these authors included fainter LRGs {i < 21) in their 
analysis that are expected to contain larger photometric 
redshift uncertainties. It is, however, not clear how the 
authors accounted for the varying redshift errors with 
galaxy brightness in their estimate. 

In contrast, we applied a 10% correction factor for 
the clustering measurements according to our simulation 
studies described in § 3.2.2. The halo masses derived 
for weak/strong Mg II are larger than the ones obtained 
bv IBouche erall (j2006l ) over similar Wr(2796) intervals. 
The apparent discrepancy may be due to the absolute 
halo bias of LRGs included in different analysis. We cal- 
culated 6hg directly from the LRG cluster ing signal ob- 
served in a volume-limited sample, whereas lBouche et al.l 
(|2006( l compared their LRG clustering strength with pre- 
vious studies that were carried out using a similar, but 
not identical galaxy population. Additional uncertainties 
in the estimated correction factor to account for photo- 
metric redshift errors may also contribute to the differ- 
ences in our findings. 

5. DISCUSSION 

We have calculated the clustering amplitude of Mgll 
absorbers with respect to three samples of LRGs. Us- 
ing the volume-limited sample, we have computed the 
absolute bias and typical halo mass of two subsamples 
of Mgll absorbers: 14^^(2796) > 1.5 A and W^(2796) = 
1 — 1.5 AA la anti-correlation between Wr{2796) and 
mean mass is seen in panel (c) of Figure|6l If a significant 
anti-correlation signal is confirmed by larger datasets, 
this would imply that weaker Mgll absorbers are found 
to be more strongly clustered than stronger ones. Our 
results show that a significant fraction of the Mgll ab- 
sorber population of Wr{2796) — 1 — 1.5 A absorbers 
are found around massive galaxies with logMh < 13.4, 
whereas absorbers of Wr{2796) > 1.5 A are primarily 
found in log Mji < 12.7 galaxies. Larger datasets for 
both LRGs and Mgll would improve the precision on 
the clustering measurements and the equivalent width 
vs. mass relationship. Corresponding galaxy luminosi- 
ties can be inferred from the bi as-luminosity relation 
found fo r SPSS data at z ^ 0.1 (jTegmark et al.l 12004 
see also IZehavi et al.l 120051 ). Assuming 6 = 1 for L,- 
galaxies (which is reasonable sinc e L<. -galaxies are f ound 
in halos of mass lO^^-^ Mq; see I Zheng et al.ll2007D . the 



bias-inferred luminosity for weak and strong absorbers 
are w 4.5^* and w 1.5L, respectively. Note that these 
values are obtained from the mean bias estimate of equa- 
tion 15 where more massive galaxies (higher luminosity) 
have higher weights (bias) than less massive (lower lumi- 
nosity) objects. These values should not be interpreted 
as the luminosity of a typical galaxy producing absorbers 
of a given strength. In addition, the relationship between 
bias and luminosity is not a one-to-one relation since 
the bias of galaxies is aflFected by the satellites within 
their dark matter halos. There is also an expected red- 
shift evolution of t he bias-luminosity relation. However, 
I Zheng et al.l (|2007| ) found that there is little evolution 
in the halo mass hosting the central galaxies between 
z ~ and z = 1. These authors found that a typical i»- 
galaxies reside in halos only a few times more massive at 
z = 1 than at z = 0. Therefore, we assume no redshift 
evolution for the bias-luminosity relation and used the 
expression found at z = 0.1 as a reasonable guess for our 
z — 0.5 sample. 

Another important aspect of the results is that the 
Mgll-LRG cross-correlation function continues to ex- 
hibit a strong signal (even after correcting for photmetric 
redshifts) down to ~ 0.3 Mpc, indicating that some 
of the Mgll absorbers and the LRGs share a common 
dark matter halo. Here we discuss the implications of 
these results. 

5.1. The bias vs 1^^(2796) relationship 

The Wr vs. mean halo mass relationship found in our 
analysis is qualitatively consistent with the previous re- 
port by Bouche et al. (2006; see also iLundgren et al.l 
[20091) . It is a 1-tT trend but it argues against the simple 
notion that more massive halos might produce stronger 
absorbers because they contain a larger volume of gas. 
Our results show that weaker absorbers are preferentially 
found in more massive halos. 

Bouche et al. (2006) measured the Mgll-LRG cross- 
correlation function using 1806 Mgll absorbers of 
WV(2796) > 0.3 A and 250,000 LRGs of i' < 21 at 
z = 0.35 - 0.8. They found that Mgll absorbers of 
Wr(2796) < 1 A appeared to be more strongly clustered 
than the Wr (2796) > 2 A ones by nearly a factor of two. 
The authors attributed the observed anti-correlation to a 
starburst outflow origin for absorbers with W,.(2796) > 2 
A, in order to explain the on-average lower halo mass of 
these absorbers. 

However, our analysis shows that absorbers with 
Wr(2796) ^1.5 are essentially unbiased with respect to 
dark matter [b = 1). This indicates that the halo pop- 
ulation probed by these absorbers is consistent with a 
random, unbiased sample of dark matter halos, and does 
not favor a specific sub-population such as starbursting 
systems. For absorbers of Wr(2796) < 1.5 A, the mean 
halo bias was found to be still higher. This large mass 
scale is at odds with previous findings that Mgll ab- 
sorbers of Wr(2796) = 0.3 — 1 A arc associated with 
L*-type galaxies (e.g., Steidel et al. 1994). 

To understand the physical mechanisms that could 
explain a Wr(2796) vs. cluste ring amplitude anti- 
correlation, iTinker fc ChCTi l2008l (hereafter TC08) de- 
veloped a halo occupation model that constrains the 
cold gas content of dark matter halos based on the ob- 
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served number density and clustering amplitudes of the 
Mgll absorbers. In the TC08 model, Mgll absorbers 
serve as a representative tracer of cool gas (T ^ 10^ K) 
in dark matter halos, and the observed anti-correlation 
arises as a result of an elevated clustering amplitude of 
Wr(2796) = 0.3 — 1.5 A absorbers due to the contri- 
butions of residual cold gas in high-mass halos {Mh > 
10^^ Mq). These massive halos are rare and are likely 
missed in small samples (e.g., Steidel et al. 1994). The 
halo occupation model represents the first empirical con- 
straint of the cold gas content across the full spectrum 
of dark matter halos and provides additional information 
for models of the growth of gaseous halos. 

5.2. Presence of cool gas in massive halos 

From Figure [31 it is worth noting the strong clustering 
of Mgll absorbers for all three LRG samples. Indeed, 
the clustering strength is comparable to the LRGs auto- 
correlation signal for the most inner bin (r^ — 0.31 h^^ 
Mpc). In physical units, this bin corresponds to 0.21 
Mpc. For the measured clustering scale, the typical mass 
scale of LRGs is roughly log(Mh/M0) « 13.2, implying 
a typical virial radius of i?vir ~ 0.35 Mpc. The virial 
radius is larger than the physical separation probed by 
the most inner bin. The strong cross-correlation signal 
is thus indicative of the presence of cool gas well inside 
the virial radius of massive galaxies. 

To further investigate the presence of cool gas in mas- 
sive halos, we examined the effects that a varying cold 
gas covering fraction (k) in LRG halos would have on 
the cross-correlation signal. To do this, we constructed 
a mock LRG catalog based on our halo occupation dis- 
tribution fits to the volume-limited sample (see bottom 
left panel of Figure 4). Using the best-fit halo occu- 
pation function, we populated the halos identified in a 
z=:0.5 output of an N-body simulation. This simulation 
is smaller in volume (400 Mpc on a side) than the 
simulation used to test the photometric redshift errors, 
in order to probe lower- mass halos. Details about this 
simulation can be found in lTinker et al.l (|2007t ). The cos- 
mology of this simulation differs from our fiducial cosmol- 
ogy (WMAPl vs WMAP5), so the large-scale bias of the 
LRG auto-correlation function is lower in comparison to 
the data, but the one-halo clustering is a good match to 
the data. 

We then simulated a mock Mgll absorber catalog 
by selecting random sightlines. Every halo of Mh > 
10^^h~^ Mq was allowed to produce a mock Mgll ab- 
sorber if the impact parameter was less than the virial 
radius, whether or not it contains a mock LRG. We then 
measured the cross-correlation between mock absorbers 
and mock LRGs for different values of k following differ- 
ent recipes. First, all halos of containing an LRG yielded 
an absorber if intersected by a sightline. Namely, all LRG 
halos have a gas covering fraction of k = 1. Then, we 
varied the covering fraction of Mgll in the mock LRG 
halos. 

The results for the two limiting cases oi k — and 
K = 1 are shown in the left panel of Figure [7| along with 
the mock LRG auto-correlation function calculated from 
the box. The right panel shows four curves correspond- 
ing to four different k values on top of our volume-limited 
cross-correlation and auto-correlation measurements for 
Wr{2796) = 1 — 5 A absorbers. The cross-correlation 



function is a probe of the relative covering fraction of 
LRG halos with respect to lower mass halos. Lower- 
ing K by the same amount for all halos, LRGs and X* 
alike, does not change the resulting cross-correlation sig- 
nal. The results from Figure 7 imply that the covering 
fraction of LRG-hosting halos must be comparable to 
that of halos that contain L^-galaxies at their centers. 

The presence of cold gas in massive halos has been a de- 
bated subject in recent numerical simulatio n studies. In 
a serie s of high-resolution SPH simulations, iKeres et al.l 
()2008[) and iBrooks et all (|2008[) examined the temper- 
ature history of gas accreted onto dark matter halos. 
They found that most of the baryonic mass is acquired 
through filamentary cold mode accretion that is never 
shock he ated to the vi r ial tem perature for halos of < 10^^ 
h-^ Mq. IKeres et all ()2008i) found that these cold flows 
are not present in massive halos typically hosting LRGs. 
At high redshift, cold flows may penetrate inside the 
virial shock of 10"'^^ M©, but this effect is highly 
redshift dependent, and is not likely to yield to cross- 
correlation functions seen in Figure 4. 

Other mechanisms such as thermal instability could 
generate pockets of cold gas ii iside a hot medium (e.g. , 
iMo fc Miralda-EscuddfT99l . iMaller fc BulT^ (pOOl 
showed that the hot gas is thermally unstable and prone 
to fragmentation. They also show that cooling proceeds 
via the formation of cold 10^ K clouds in pressure equilib- 
rium with the hot halo gas. For a Milky- Way-size system, 
cool clouds of mass ~ 5 x 10^ Mq are expected to extend 
up to ~ 150 kpc from the galactic cente r and survive 
for several Gyrs. iKaufmann et al.l (j2008[ ) showed that 
cloud formation is viable in halos, but needs to be ex- 
te nded to higher mass. In a sim ilar argument developed 
bv lMo fc Miralda-Escudd (|1996( ). a two-phase medium in 
pressure equilibrium was used to explain observations of 
Lyman limit systems. This model also makes predictions 
about the presence of C IV around low-mass galaxies and 
at large impact parameters of massive gala xies. These 

JTedictions were partially confirmed later by IChen et al.l 
2001h who found C IV in galaxies of different morpholo- 
gies and luminosities. These authors also observed the 
sharp boundary in the W r-(pro,iected separation ) plane 
that was also predicted. iMo fc Miralda-Escudd ()1996l ) 
attributed the presence of Mg II absorbers to cold pock- 
ets of photoionized gas in halos around massive galaxies. 
The observed strong Mg II-LRGs cross-correlation signal 
on scales smaller than the virial radii of halos hosting 
LRGs indicates the presence of cold gas is more common 
around massive galaxies than previously thought. 

On the observational side, the detections of cool gas in 
group size halos is uncertain. The challenges lie in the 
limited sensitiviti es available to detect HI gas via 21-cm 
observations fe.g..lVerdes-Montenegro et al]|2001l . [20071) . 
IVerdes-Montenegro et al.l ( 2007| ) reported detections of 



HI column density down t o ~ 10 cm~ in some of the 
Hickson Compact Groups (|Hicksonlll982[ ). Observations 
are not sensitive enough to probe the low HI column 
density environment yet. 

The difficulty in finding cool gas around groups using 
21-cm observations underscores the powerful application 
of QSO absorption- line studies. We are currently con- 
ducting a follow-up imaging and spectroscopy campaign 
to study the cool gas content of individual LRG halos 
(Gauthier et al. 2009 in preparation). 
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Fig. 7. — Left panel: Cross- and auto-correlation results for the 400 h ^ Mpc box. The solid line shows the auto-correlation of the mock 
LRG sample. The diamond points show the cross-correlation results for which k = 1 for all halos above lO^'^ Mq. The square points 
represent the case of no cold gas in the halos of the mock LRGs and k = 1 for all remaining halos > 10^^ Mq. Right panel: the cross- 
and auto-correlation results for the volume-limited sample of LRGs (Wr(2796) = 1 — 5 A) are shown in gray and black points respectively. 
The four curves correspond to different values of k. for the mock LRGs. From bottom to top, ft = 0,0.33,0.5, 1. 



5.3. Future prospects 



DR7 Mgll database & HOD 
modeling 



Figure 7 demonstrates how a halo occupation approach 
to the galaxy-absorber cross-correlation function at both 
large and small scales can put constraints on the cov- 
ering fraction of cold gas in LRG-hosting dark matter 
halos. In a forthcoming paper, we will address the de- 
tailed halo occupation distribution modeling of the Mg II 
absorber environment with a particular focus on the one- 
halo term. This will be achieved by using the SDSS DR7 
Mgll absorber catalog. This catalog, currently in prepa- 
ration (Prochaska et al. 2009), will more than double the 
number of absorbers from the current DR5 sample. It 
will allow us to probe smaller projected separations and 
improve the clustering measurements at rj, < 1 Mpc. 
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